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Abstract 


In the solid desiccant wheel air-conditioning system, the performance of the desiccant wheel is 
critical to the capability, size and cost of the whole system. Constructing mathematical model is an 
effective method for analyzing the performance of the rotary wheel as well as the system. The model 
can also be used to guide system operation, interpret experimental results and assist in system design 
and optimization. The overall objective of this paper is to provide a review of various efforts that 
researchers have made to mathematically model the coupled heat and mass transfer process 
occurring within the wheel. The paper first briefly describes desiccant wheel including fundamental 
principle, heat and mass transfer mechanism and the method of model establishment. Then various 
models consisting of ideal assumptions, governing equations, auxiliary conditions, solution methods 
and main results are presented. The models can be classified into two main categories: (1) gas-side 
resistance (GSR) model; (2) gas and solid-side resistance (GSSR) model which can be further 
subdivided into pseudo-gas-side (PGS) model, gas and solid-side (GSS) model and parabolic 
concentration profile (PCP) model. It shows that GSSR models are higher in precision and more 
complex compared with GSR models. In addition, the simplified empirical models based on 
measured data are briefly discussed. This review is useful for understanding the evolution process and 
status quo of the mathematical model and highlighting the key aspects of model improvement such 
as taking account of pressure loss or air leakage. 
© 2007 Elsevier Ltd. All rights reserved. 
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1. Introduction 


Most of traditional vapor compression air conditioning system influences the 
environment in three aspects. First, the CFCs and HCFCs used as working fluids are 
ozone-depleting substances. Second, most of the refrigerants such as HCFCs are the global 
warming gases. Third, refrigeration system consumes energy which in general is produced 
from fossil fuel and related to release of pollutions. These polluting materials such as CO, 
SO,, NO, are harmful to the health of people. While desiccant wheel air conditioning 
systems offer a lot of significant advantages over other AC systems. A few remarkable 
advantages can be summarized as: 


1. Water, a natural working fluid, is used as refrigerant. Desiccant materials such as silica- 
gel and zeolite are also environmental friendly. 

2. Desiccant cooling system can be powered by low grade thermal energy sources such as 
solar energy, geothermic energy and waste heat. Hence the operating costs can be 
significantly reduced. 

3. The system deals with the sensible and the latent loads, respectively, and can be utilized 
in various climates or environmental conditions. 

4. It overcomes the discontinuous problem of the fixed bed desiccant cooling system. 

5. Compared with the liquid desiccant system in which the liquid and air directly interact, 
the solid one is compact and less subject to corrosion. 
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Nomenclature 


Do 
Ds 

f 
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BNET tN 
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the distance in radial direction (m), Eq. (20) 

channel cross section area (m°) 

transfer area per unit volume (m*/m*), Eq. (65) 

outside atmospheric pressure (Pa), Eq. (26) 

specific heat (J/kg K) 

hydraulic diameter de = 44/P, Eq. (83) 

desiccant effective diffusivity (m7/s), Eq. (123) 

gas phase diffusivity (m7/s), Eq. (116) 

Knudsen diffusivity (m/s), Eq. (132) 

ordinary diffusivity (m/s), Eq. (132) 

surface diffusivity (m*/s), Eqs. (108), (131) 

desiccant fraction in the wall or in the matrix 

the mass of support materials/desiccant per length in a tube (kg/m), Eq. (32) 
ratio of free flow area to cross section area of rotary wheel, Eq. (65) 
air-side convective heat transfer coefficient (W/m? K) 
compound convective heat transfer coefficient (W/m? K), Eq. (83) 
convective mass transfer coefficient (m/s), Eq. (47) 

heat of vaporization of water (J/kg) 

enthalpy (J/kg) 

thermal conductivity (W/m K) 

gas-side mass transfer coefficient (kg/m? s) 

composite mass transfer coefficient (kg/m? s), Eq. (84) 

actual process channel length (wheel thickness) (m) 

Lewis number 

mass (kg) 

mass flow rate of the stream (kg/s) 

rate of phase change per unit length (kg/ms), Eq. (42) 

mass flow rate of air per unit cross section area (kg/m7s), Eq. (65) 
the number of the channel, Eq. (20) 

number of transfer units, Eq. (91) 

channel cross section perimeter (m) 

heat of sorption (J/kg adsorbate) 

radial displace or radius of the particle (m), Eq. (143) 

radius of the matrix axis (m), Eq. (100) 

radius of the rotary wheel (m), Eq. (100) 

separation factor (reflecting desiccant isotherm shape), Eq. (16) 
temperature (K) 

time (s) 

time required for a complete matrix rotation (s), Eq. (82) 

air velocity (m/s) 

the total volume of desiccant wheel (m°), Eq. (107) 

desiccant adsorption mass (kg adsorbate/kg adsorbent) 

actual channel width (m), Eq. (20) 

humidity ratio (kg/kg) 
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absolute humidity ratio of the air stream (kg water vapor/kg dry air) 
humidity ratio in equilibrium with the desiccant (kg water vapor/kg dry air) 
absolute humidity ratio of the air stream (kg water vapor/kg moist air) 
humidity ratio in equilibrium with the matrix (kg water vapor/kg dry air) 
axial displace through matrix measured from period entrance (m) 


Greek letters 


p relative humidity 

n fraction of phase change energy delivered directly to the air, Eq. (42) 
p density (kg/m°) 

o volume fraction, Eq. (53) 

ô the thickness of desiccant felt (m) 

0 time from the beginning of period (s), Eq. (82) 
0; duration of period j (s), Eq. (82) 

(0) wheel rotate speed (r/h) 

T tortuosity, Eq. (109) 

Er felt porosity, Eq. (131) 

My dry matrix mass divided by dry fluid mass in the matrix in period j, Eq. (79) 
Subscript 

0 initial estate 

1 the process air 

2 the regeneration air 

a dry air 

ad air in desiccant pore, Eq. (116) 

d desiccant 

eff effective property, Eq. (140) 

f entire felt (matrix and pores), Eq. (139) 

fd dry solid felt material, Eq. (138) 

fm matrix part of felt, Eq. (138) 

Ip pores (inter-particle voids), Eq. (138) 

fs dry solid felt material, Eq. (138) 

g gas including dry air and vapor 

in inlet condition, Eq. (5) 

j period index (during the period j), Eq. (77) 

l liquid water of frost, Eq. (52) 

le local, Eq. (37) 

m matrix 

s surface, in gas phase adjacent to gel particles, Eq. (130) 
sp support materials, Eq. (30) 

v water vapor 

vs saturated water vapor, Eq. (26) 

w wall, Eq. (119) 

ws saturated water, Eq. (75) 
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Desiccant wheel found in the desiccant cooling system is widely acknowledged as the most 
critical part of the system, both in terms of cycle performance and system cost. It has been 
successfully used in many systems for dehumidification and enthalpy recovery. With the 
development of computer hardware and numerical methodology, mathematical models are 
being used to carry out critical investigations concerning on the desiccant wheel. The 
advantages of this method are that it takes less time and cost than experimental method for 
predicting the performance of a desiccant wheel. The second point is that the mathematical 
model can produce large volumes of results at virtually no added expense and it is very 
convenient to perform parametric research and optimization analysis. The third reason is 
that by experimental means, some parameters like that inside the tubes cannot be 
measured. The last but not the least reason is the fundamental physics and surface 
chemistry of rotary desiccant wheel, fixed desiccant equipment and enthalpy wheel are 
similar, so that other installations may be investigated too by utilizing the model of 
desiccant wheel [1-3]. Consequently, constructing valid mathematical model for desiccant 
wheel has become the key subject of many studies. Lots of mathematical models have been 
constructed and employed to analyze, develop and design desiccant wheels. However, to 
the authors’ knowledge, there is not yet a work that summarizes comprehensively the 
numerous significant works done on modeling the desiccant wheel system. Hence, the 
objective of this article is to examine the progress in the area of mathematical modeling on 
the desiccant wheel. 

Based on the heat and mass conservation principles, the combined heat and mass 
transfer in desiccant wheel can be formulated with complex mathematical detail as 
governing equations and auxiliary relations. This is the main and universal method for 
constructing mathematical model. There are also some simplified models for ease of 
computation such as pseudo-steady-state model [4,5], which are not fundamentally correct 
or universal. Moreover, it is an effective way to get the empirical performance fitted 
correlation using large amounts of experimental results. But its application is restricted in 
the specific range of experimental working conditions. The review of this paper will focus 
on the mainstream mathematical models. Besides, some empirical models will be briefly 
introduced too. 

The first category of the paper is distinguished by whether the model gives solid-side 
resistance (SSR). Within the first section, it is ordered according to the dimension of 
models. For the second category, papers are grouped together by the method of calculating 
SSR. For each model reviewed, attempt has been made to gather up the assumptions, 
governing equations, auxiliary relations and the solution techniques for ensuring the 
integrality of the model. The results and validation of some models are listed when 
comparison between them are needed. 


2. Problem description 
2.1. The fundamental principle 


Desiccant wheel is also known as dehumidification wheel, wheel dehumidifier or rotary 
dehumidifier in different literatures. The main components of desiccant wheel include 
matrix, clapboard, the wheel case, the air heater, the driven motor, etc. as shown in Fig. 1. 
Matrix usually contains supporting material and desiccant material. The case of desiccant 
wheel is driven by an engine and moves at a given rotary velocity. It consists of a large 
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nonhygroscopic 
and adiabatic ektexine 


Felt/Supporting material 


Adsorbent 
one channel 


Fig. 1. Schematic figure of rotary desiccant wheel. 


number of channels whose walls are constituted by supporting material and coated or 
impregnated with desiccant material. The cross section of wheel is divided into process air 
side and regeneration air side by the clapboard. Water vapor is adsorbed by the desiccant 
when moist air is passing through process air side. And while the hot air which is heated up 
by the heater flows through regeneration side, water is desorbed out from the desiccant and 
desiccant material is regenerated. The airflows can be either in counter-flow or co-current 
arrangement. It can be found that the reduction of moisture in the process air is at the price 
of the energy consumed to heat the regeneration air. This mode is called active desiccant 
wheel. If no air heater is built in system, the mode is named as passive desiccant wheel, 
enthalpy wheel or rotary energy wheel. It should be noted that the wheel is usually installed 
with thermal insulation and air proof material. Hence the energy and mass exchange 
between the wheel and environment is neglected in most cases. 


2.2. Mechanism of heat and mass transfer in desiccant wheel 


During adsorption, heat of sorption is released as desiccant adsorbs water vapor. The 
heat generated in the desiccant is transmitted through the material which decreases the 
sorption capacity. Therefore, the heat and mass transfer within the desiccant are coupled 
and should be considered simultaneously in developing mathematical model. 

In a desiccant dehumidifier, there exist gas-side resistances (GSRs) and solid-side 
resistances (SSRs). GSRs come from diffusion and heat conduction within the air and 
convective heat and mass transfer between air and desiccant. It influences the heat 
and mass transfer from the air stream to the surface of the desiccant felt. Due to diffusion 
and heat conduction within the air is small in comparison with convective transfer, they are 
always neglected. SSRs come from heat conduction and mass diffusion within solid 
desiccant felt. Where the mass diffusion involves a number of physical mechanisms: (1) 
surface reaction combined with the adsorption process; (2) interparticle channel diffusion 
and adsorption within the layer; (3) micropore and macropore diffusion and adsorption 
within the desiccant particles [3]. 
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2.3. Mathematical model establishment 


The main steps to establish mathematical model include: (i) assumption proposal and 
selecting appropriate control volume; (ii) deriving the partial differential governing 
equations based on the mass and energy balances; (iii) providing auxiliary relation to close 
the governing equations; (iv) adopting proper mathematic approach to solve the model. 

Some ideal assumptions are needed for the reason that: 


(1) The mass and heat transfer occurring in the wheel are too complicated to understand 
completely. 

(2) Some factors that do not influence heat and mass transfer significantly can be 
neglected. Then the complexity of the governing equations can be reduced and the 
solution is relatively simple and time-saving. 


2.3.1. General governing equations 

Euler cylindrical coordinate system has three coordinates: axial direction (z), radial 
direction (r) and circumferential direction (@) and is usually employed to describe the 
wheel which has the same cylindrical geometry like Euler cylindrical coordinate system. In 
general, governing equations involve four main equations such as: mass and energy 
balance equations in the air and in the solid side. For clarity, the system is considered to be 
one-dimensional and explanations of main terms in governing equations are presented 
below. The physical model is shown in Fig. 2. 

The moisture conservation in the air can be expressed as 


Ot Oz 


The first term on the left-hand side (LHS) of the above equation is the moisture storage 

term in the air. The second term on the LHS represents the rate of moisture variation in the 

air due to the axial flow of the air. The first term on the right-hand side (RHS) of Eq. (1) 

expresses the rate of moisture variation in the air caused by the convective mass transfer. 
The energy conservation for the air is described as 


OY, OY, 
depa (= +u ) = K,( Ya- Ya). (1) 


OTa oTa ka a? Ta 
dacroa( ar +u Ag CoaPa zZ ) = h(Ta E Ta) F Cov Ky Ya- Yal(Ta = Ta). 


(2) 


Process air 


et ie air 
Moist air 


Fig. 2. Euler cylindrical coordinate system and a model of control volume of the desiccant wheel. 
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The first term on the LHS of the equation represents energy storage in the humid air. The 
second term on the LHS calculates the rate of energy variation in the air due to the axial 
flow of the air. The third term on the LHS represents the heat conduction in the air. The 
first term on RHS of Eq. (2) expresses the convective heat transfer between the air and 
solid desiccant. The second term on RHS of Eq. (2) indicates the sensible heat transfer 
between the air and solid desiccant. 

The moisture conservation in desiccant is given by 


De 


Ow W 
Pao\ y 622 


) = K,(Ya— Yu): (3) 


The first term on the LHS of the above equation is the moisture storage term inside 
desiccant. The second term on the LHS shows mass diffusion within the solid desiccant in 
axial direction. The first term on the RHS of Eq. (3) represents the convective mass 
transfer between the air and desiccant. 

The energy conservation in desiccant can be written as follows: 


OTe ke OTe 
cra8( at CpyPa 072 — WT = Ta) = Ky( Yq ~ Yast 


F CwKy(Ya = YaX(Ta = Ta). (4) 


The first term on the LHS of the above equation is the energy storage term of desiccant. 
The second shows the heat transfer due to heat conduction within desiccant. The first term 
on the RHS of Eq. (4) calculates the convective heat transfer between the solid desiccant 
and air. The second term on the RHS expresses the influence of adsorption heat. The third 
term on the RHS expresses the sensible heat transfer between the solid desiccant and air. 

It is important to note that in the above governing equations, it is assumed that all of the 
adsorption heat is delivered to desiccant. 


2.3.2. Auxiliary conditions and mathematical methods 

In order to close the partial difference equations, some auxiliary conditions have to be 
provided. They are boundary conditions, initial conditions for transient analysis, 
adsorption equilibrium correlations. Also, the values of adsorption heat qss convective 
transfer coefficient h, thermal conductivity k and moisture diffusivity D, in desiccant 
material are need to give. 

a. Boundary conditions: Boundary conditions describe the behavior of the simulation at 
the edges of the simulation region. The boundary conditions, commonly encountered in 
desiccant wheel model, are usually in three modes: 


(1) Boundary conditions specify the values or the normal derivative of parameters on a 
surface [6], such as: 
In the process/regeneration air (i = 1 for process air; i = 2 for regeneration air) 


Tain = Tijin, (5) 


Yain = Yin. (6) 
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The insulated and impermeable boundaries result in [7] 


0 Tn, 0 Tm 0 Tm 
= = =0, (7) 
Or |o OZ | 25 Oz |r 
0 Yn ð Yn 0 Yn 
= = = 0. (8) 
or r=0 Oz z=0 0z z=L 


(II) Periodic boundary conditions: The periodic boundary conditions at ® = 0 [8] 
Y,(0, 2, t) = Y,(27, z, t), 
T (0, z, t) = T27, z, t), 
W(0,z, t) = W(2z, z, t), t 
Ta(0, z, t) = Ta(2r, z, t). 
The periodic equilibrium boundary conditions expressed in the form of limitation 


[9]: 


for 0<z<L 


ae [Wagn] = lim, [Wn(1-F0)], 
jon [Hn (5 n) = Past [Hm (1 = T n)| > 
l z to Ik z (10) 
i, [wal] = hy ef) 
i, ta) =n [P(E] 
(II) Boundary conditions depict heat and mass transfer on a surface [10] 
OYa ow 
PppDe az (2,0) + paDs = (62,0) = 0, (11) 
OYa Ow 
-ppDe A(t, 2,44) — paDs y (t204) = K[Y a(t, zaa) — Yalt2), (1D) 
OTa oT, 
ka—(t = Py OCpw ——» 1 
d or (t, z, 0) Py Cp Ot ( 3) 
oT 
ka ZA (62, aa) = WT a(t, z, aa) — Talt, 2). (14) 


b. Initial conditions: Initial conditions give the value of parameters at a given starting 
time (t = 0). For example, the initial condition of parameter X is 


X(t = 0) = Xp. (15) 


X is replaced by the detailed parameters in different models. 
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c. Equilibrium adsorption equation: The equilibrium equation determines the amount of 
moisture that will be adsorbed by the desiccant. Emphasis should be placed on the fact that 
adsorption is a dynamic process. With different operation conditions, such as relative 
humidity of the air in equilibrium with the adsorbent ®, desiccant temperature Ty, etc., 
desiccant has different adsorption capacity and equilibrium adsorption quantity. Up to 
now, through experiments to obtain the dynamic adsorption data and equilibrium 
equation for different desiccants is the main method. As a result, there are many types of 
equilibrium adsorption relationships including: relative humidity of moist air in 
equilibrium with the adsorbent ® and separation factor R are used to calculate the water 
content of the adsorbent W: 


W/Wmax = D/[R + (1 — R)9]. (16) 


Some researchers consider that W is the function of ® and temperature of moist air in 
equilibrium with the adsorbent Tg: 


W = f(Ta, Ya) =f (Ta, ®). (17) 


Sometimes it is assumed that the Wis a function of ® only and invariable with temperature 
Ti 


W = f(®). (18) 
Besides, the following relationship is utilized in many works: 
® = f(W, W max, R, h*, Pys). (19) 


d. Mathematical solution: Due to the complexity and nonlinearity of such coupled heat 
and mass transfer models, it is difficult to obtain analytical solutions. Hence numerical 
methods have been applied to solve the partial differential governing equations. At 
present, the most widely used solution method for desiccant wheel model is finite difference 
method. Other methods such as finite volume method [7] and analogy method [11,12] have 
also been used in some works. 


3. Gas-side resistance (GSR) model 


In early research, due to the lack of study on complex diffusion mechanism inside solid 
desiccant, the heat and mass transfer processes from air stream to desiccant are thought to 
be controlled by GSRs only. Gradients of temperature and concentration within solid 
desiccant are ignored. Thus the second terms on the LHS of Eqs. (3) and (4) are omitted. 
Some researchers only consider the transfer processes in axial direction. Others consider 
that in axial as well as circumferential direction. 


3.1. One-dimensional model 


To save the computation time, it is desired to take into account as less dimension as 
possible in the model. Since the variation in moisture as well as temperature is more 
significant in axial direction, it is often chosen as main direction for heat and mass transfer 
in constructing the governing equations. Usually, condensation is not considered in the 
models. However, it happens in certain conditions. Hence some researchers include it in 
their models. 
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3.1.1. Models without considering condensation 
Charoensupaya and Worek [13] have derived a simple one-dimensional dehumidifier 
model based on the following ideal assumptions: 


1. The axial direction (z) heat conduction and molecular diffusion within the desiccant 
coating were negligible. 
2. There were no radial direction (r) temperature or moisture content gradients in the 
matrix. 
3. The axial direction heat conduction and water vapor molecular diffusion in the air 
were negligible. 
. The inlet air conditions were uniform radially, but might vary with time. 
. The tubes that make up the rotary desiccant wheel were identical. 
. All the channels were assumed to be adiabatic and impermeable. 
. The effect of centrifugal force was neglected. 
. The thermophysical properties of the dry air and the properties of the dry desiccant 
material as well as the matrix were constant. 
9. The pressure and velocity losses of the air stream in the axial direction were negligible. 
10. The heat and mass transfer coefficients between the air stream and the desiccant wall 
were constant over the entire dehumidifier. 


ooann A 


It can be inferred from the assumptions 1 and 2 that this model neglects the SSRs. Because 
of symmetry characteristic, the middle of one tube could be considered adiabatic and half 
size of the channel was selected as the differential control volume as can be seen in Fig. 3. 

The channel was assumed to be adiabatic and impermeable, hence the conservation of 
moisture and energy between the process air stream and matrix were 


oY Ma OY Ow 
{ = 4 2ô—— = 0, 20 
Pala ay + ox Oe + F0m20 > (20) 
mng (10H, . OH, My CHm 
: =0; 21 
Xin G Ot Oz XmL Ot 2) 
Conservation of moisture in the matrix was 
fmm OW 
—— =2K,(Y,—- Ym), 22 
NXmL Ot oh ) 2) 
sorbent desiccant Paan & impermeable wall 
26 
Adsorption air Desorption air 
—> —;— -mm mm MS "a, NE 
! T, Y, i 
l Y 
T, W 


Fig. 3. A model of dehumidifier differential element for GSR model [13]. 
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where m,, is the total mass of the matrix in the dehumidifier 
Mn = 26p,,nX pLa (23) 


Conservation of energy in the air was calculated by: 


mg (10H,  OHy oH, 
= 2K, Ym— Ya 2h Tm — Ta), 24 
NX m G Ot a Oz vl Joy, T a ) 4) 
where 0H,,,/OW means the differential heat of wetting and 0H,/OY, is defined as 
0H, 10H», 
= = 2 
OY, TF ow en 


The subscript “m” in the above equations means the matrix including nonhygroscopic 
exterior surface, support materials and the desiccant. Parameters n, X„ and L represent 
holistic configuration of the wheel. The governing equations were solved subject to the 
equilibrium relation, Eq. (16). It was supposed that the variation of adsorption heat q,, did 
not influence system performance greatly, hence q,, was taken as constant and the average 
value of 2700 kJ/kg HO was used. 

Then the models were converted into nondimensional form and solved numerically using 
an explicit finite-difference method. This model had been used to calculate an open-cycle, 
adiabatic, solid desiccant cooling system operating in the ventilation mode. The system 
used an adiabatic dehumidifier, a regenerative heat exchanger whose effectiveness was 
0.95, and two direct evaporative coolers with effectiveness of 0.95. The effects of 
dehumidifier channel length, desiccant mass fraction and desiccant isotherm shape on 
system performance in terms of COP were investigated. The results showed that there are 
optimum mass fraction and desiccant isotherm shape which will result in an optimum 
COP. Furthermore, the desiccant cooling system operating in the ventilation mode can 
yield COP as high as 1.4. 

Zheng and Worek [14] have established a transient one-dimensional desiccant wheel 
model which has similar assumptions and control volume as those of Charoensupaya and 
Worek [13], with certain modifications. They omitted the storage terms 0Y,/0t and 0H, /0t 
in Eqs. (20), (21), (24) and adopted the enthalpy of water vapor to replace 0H,/0Y, in 
Eq. (24). 

The new governing equations were subject to boundary condition of Eqs. (5), (6) and 
initial conditions with parameter Y,, and 7,,. The relationship between humidity ratio and 
relative humidity was described as 


0.62198P, _ 0.62198b, 


= = 26 
4 B-P, B/Pys — Ga’ 6) 
where P,,, presents the saturation pressure of water vapor 
Pws = exp Ti Ay + A3T +44T? + 15T? + ma 
T ° ln TJ’ 
A, = — 5800.2206; Az = 1.3914933; 23 = —0.04860239, a 
4 = 0.41764768 x 1074; As = —0.14452093 x 1077; 46 = 6.5459673. (0 
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The following empirical form was used to give heat of sorption [15]: 


exp W / W max) — exp(A) 


a = h*(W/Wmax) = 1 + Ah E 2, = constant, (28) 
where hy is the heat of vaporization of water. 
The equilibrium equation was represented by 
WR PD) 
p= 29 
W max F (R a D)W eo C ) 
And the enthalpy of desiccant bed was expressed by 
w 
Hm= (d — f )CpspTm +floutn+ f (H= a) aW, (30) 
0 


The simulation aimed to study different numerical methods. They adopted an implicit 
scheme finite difference method with variable step sizes in time and space coordinates to 
solve the governing equations. It has been reported that compared with explicit scheme, 
the implicit scheme is unconditionally stable, is of high-order accuracy, and enables more 
rapid calculation of dehumidifier performance with much larger time and space steps. 
Meanwhile, the simulation program was validated by two methods: (i) varying step sizes in 
both space and time; (ii) comparing with other available works, such as results of Ref. [13]. 
Both validation methods verified the convergence and accuracy of the model. Using this 
model, the effect of rotation speed on the performance of the adiabatic desiccant wheel was 
parametrically studied, and the optimal rotational speed was discovered by examining 
outlet humidity profiles of the adsorption-side as well as humidity wave fronts inside the 
desiccant dehumidifier. 

Moreover the same model was used: (i) to optimize the wheel rotational speed [16]; (11) to 
investigate the influences of desiccant sorption properties, heat and mass transfer 
characteristic and the size of the wheel on the dehumidification performance of the 
desiccant wheel [17]; (iii) and to investigate the wheel performance at the real operating 
conditions which vary significantly from the design point (i.e., at ARI conditions) [18]. 

Zhang et al. [6,19,20] constructed a mathematical model of a honeycombed desiccant 
wheel. The tube with honeycombed geometry and differential control volume are shown in 
Figs. 4 and 5. Heat of sorption was considered to be delivered to the matrix only. And the 
channel was assumed to be impermeable but not adiabatic. Hence, the model consisted of 
conservation equations of moisture and heat in the air (Eqs. (31) and (33)) and within 
desiccant (Eqs. (32) and (34)) separately, which are different form the governing equations 
(20) and (21). 


OYa OY, 
ao +u at) =K, Pu- Yo) (31) 


1, ow 


Sia = K,P(Ya— Ya), (32) 
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Air duct 


: ; Computation domain 
Matrix material and 


desiccant layer ~ \ 


Fig. 4. Corrugated air duct and the computation domain of a GSR model [19]. 


Solid desiccant 


Fig. 5. The differential control volume in the corrugated duct [19]. 


OT, OTa 
APalCpa F taal ar +u az ) = hP(Ta = Ta) F K,Pcp( Ya- YalTa = Ta), 
(33) 
1 OTa 
zV alpa + Wepi) + f sppsp] a 
= AP(Ta — Ta) + KyP(Va — Yas: + KyPepo(Va — YaXTa — Ta). (34) 
The adsorption heat of RD silica gel was calculated by the following equation [15]: 
Ast = hy [1.0 + 0.2843 exp(— 10.28 W)). (35) 


The equilibrium relative humidity on the surface of RD silica gel was obtained by fitting 
fourth-degree polynomials to the test data provided by manufacturer [21] 


® = 0.0078 — 0.05759 W + 24.16554W? — 124.478 W° + 204.226W". (36) 


The relationship between humidity ratio and relative humidity was described as Eq. (26). 
The pressure of saturated water vapor P,,, was given by Antonine equation: 


3816.44 ) 


Ta — 46.13 GY) 


Pys = exp (23.196 — 
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In addition, convective heat and mass transfer coefficients were determined by [22] 


ckaP 

h= L (38) 
Sh DP 

Ky = pa > (39) 

where 
Nuc = Shr, (40) 
P T 1.81 
D, = 2.302 x 10 (=) , Py =0.98 x 1075P,, To = 256K, (41) 
a 0 


where D, is the mass diffusion diffusivity of water vapor in the air. 

The model was solved by backward implicit difference scheme and validated by 
experimental results of a real honeycombed paper-based silica gel desiccant wheel. 
A good agreement was obtained between the calculation results and the experimental 
data. Using the model, the temperature and humidity ratio profiles in the wheel during 
both the dehumidification and regeneration processes were analyzed [20]. The results 
indicated that a hump curve of air humidity ratio along the channel existed all the 
time in the regeneration side. The performance of desiccant wheel was significantly 
affected by the hump curve. In order to improve the performance of the desiccant 
wheel, it was essential to accelerate the hump moving from the duct entrance to the duct 
exit as quickly as possible. Also the effects of velocity and inlet temperature of regeneration 
air and velocity of process air on the hump moving speed were investigated too. It was 
found that increasing the velocity and inlet temperature of regeneration air could 
accelerate the hump moving. Meanwhile the velocity of process air had little influence on 
the hump. 

Simonson and Besant [7] have presented a numerical model for coupled heat and 
moisture transfer during adsorption and desorption processes. Although proposed for 
rotary energy wheel, the model is still listed here because: (1) desiccant wheel and rotary 
wheel have similar operation and function of transferring heat and moisture; (2) a new 
factor y is proposed in their model to account for the distribution of phase change energy 
between the desiccant and the air. The tube geometry cross-section and limited simplifying 
assumptions are basically similar to those of Zhang [6]. 

One tube is chosen as the control volume (Fig. 6). The conservation of energy in the air 
is given by 


aT, ar, 
PglpaA + UPylygA ae = gyn +hP(Tm —T,), (42) 


where y is the fraction of the phase change energy entering the air. And then the 
first term on the RHS of above equation calculates the heat of sorption delivered to the air 
side. 
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we & adiabatic surface 


desiccant 


r air 


Z 
=P (axial coordinate) 


Fig. 6. Side view of rotary energy exchanger duct showing the aluminum surface [7]. 


The conservation of energy in the matrix (subscript “m”, including aluminum, desiccant 
and moisture) is described by 


OTm ô OT m . af 
PinCpmAm -57 = az (KA a) = m1 = n) +m cpi(Tg ~ Tm) F hP(T, ~ Tm). 


(43) 


It is supposed that the heat conduction through the matrix in the axial direction is 
dominated by the conduction through aluminum surface layer (subscript “al”) which is 
expressed by the second term on the LHS of Eq. (43). The first term on the RHS calculates 
the rest heat of sorption delivered to the matrix. The second on the RHS shows the sensible 
heat transfer between the solid desiccant and air. 

Conservation of mass in the air results in two equations, one for the dry air and another 
for water vapor, which is different from the general governing equation (1): 


õp, , 9 ss 
A art By (Pot) =M, (44) 
op 0 
a4 — = 0. 4 
Pa + = (pH) =0 (45) 
Conservation of mass in the matrix is given as 
ow, 
Pdary4d ET m, (46) 
where “m” means the rate of phase change per unit length 
i! = hmP(p, ~~ Poy)» (47) 


where p, and p,m stand for water vapor density and water vapor density on the desiccant 
surface. The following equation is used to calculate the value of py m: 


Pis Ti W C 
— Pus(Tm) 0<W<Wn——, 
ENS 


; 4 
Pon R,T m (W m = W)C’ ( 8) 
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where P,, is saturation vapor pressure and calculated by 


Pos = exp(F), 

C 

a C+ GT + CGT? ter + CTt+ Giant 173K<T<273K, 
Fata 

FT Co + CoT + CT? CRT + CplnT, 273 K<T <473 K. 

(49) 
Heat of sorption is given by 
W A 

Ist L h*(W JW max) = 1 + w(i = ) . 2=f0), (50) 
hy W max 


where 4 describes the shape of heat of sorption as a function of moisture content. 
The geometry and property relations for the moist air (g), desiccant (d) and matrix (m) 
are 


CpgPg = Palpa + Py pv» (51) 
Pa = Paay W +1), CpdPa = Paary(Wepi + Cpary), (52) 
Pm = FdPg + Capas CpmPm = Pm = FdPalpd + FalPaiCpal, (53) 
Ca + 0a = 1. (54) 


The convective heat transfer coefficient A in Eqs. (42), (43) was obtained from numerical 
simulations of simultaneously (thermal and hydrodynamic) developing flow in an 
equilateral triangular duct [23]. Convective mass transfer coefficient A„ in Eq. (47) was 
determined by using the analogy between heat and mass transfer. The boundary conditions 
for the problem included Eqs. (5) and (7) and 

Pain = Piin Poin = Pipins Yin = Uiin- (55) 

The governing equations were discretized by finite volume method with a staggered grid 
[24]. Upwind difference scheme and central differencing were used in the air and matrix, 
respectively. Transient terms and source terms were treated in fully implicit and semi- 
implicit manner. The discretized equations were solved using a Gauss-Seidel iteration 
technique with under relaxation. 

The model was validated by comparing its results with the experimental data of a 
commercial molecular sieve rotary energy exchanger by Simonson and Besant [25]. It has 
been reported that the measured and simulated results agreed well within experimental 
uncertainty. Besides, the effects of some important factors, such as phase change term q, 
adsorption heat and isotherm, axial heat conduction in matrix, storage in air and entry 
length, on the predicted performance of molecular sieve energy wheel were discussed. The 
predicted effectiveness was found to be greatly influenced by the value of n. 


3.1.2. Models considering condensation 
In the above models, the desiccant wheels are formulated at or near atmospheric 
pressure, where no condensation happens. When desiccant wheel works under high 
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pressure, condensation might occur in the regeneration process. The models reviewed in 
this part are capable of properly simulating desiccant wheel operating at high pressures. 

In order to analyze the performance of an adiabatic desiccant wheel under high process 
and regeneration air stream pressures, Mihajlo and William [26] have derived a model 
containing two groups of governing equations, one for the noncondensation situation and 
another for the case when condensation occurred. The ideal assumptions, governing 
equations for noncondensation situation and auxiliary relations are similar to that of 
Zheng and Worek [14]. A new set of governing equations were constructed when 
condensation happened such as: 

Conservation of moisture in the process air: 


OY, 2h(T, — Tm) 
Palp -z = 2Ky(Yn — Ya) 7, i (56) 
á ‘fg 
Conservation of moisture in the matrix: 
0 W 2h Ta — Tm 
F Pm Agr S 2K, ( Ya Yn) + ( ) š (57) 
ot hy 
The temperature of air which was saturated with water vapor was taken as constant 
Ta = constant. (58) 
Conservation of energy in the matrix 
OH m 2n(T 4 = Tm) 
Pmô -Agp = 245K y( Y, = Ym) ag ltz + (59) 
Ot hy 


An implicit finite-difference scheme was adopted to solve the model. The performance 
analysis using this model revealed that: (1) the value of optimal separation factor R under 
high pressure is greater than that under atmospheric pressure; (2) if the regeneration 
temperature T,eg is between 65 and 100°C and R is between 0.01 and 0.1, condensation 
could occupy 21—40% of regeneration section of the wheel; (3) only for R larger than 0.4 
does condensation not happen. 

A similar extension of the model [6] was conducted by Zhong [27], to include 
condensation process. Similarly, two sets of governing equations were established, 
respectively, for the condensation process and noncondensation process. The assumptions, 
governing equations of noncondensation status, auxiliary relations and solution method 
were similar to that of Zhang [6]. The derived model of condensation process consisted of: 

Conservation of moisture in the process air: 


OY. OYa\ _ hP(Ta — Ta) 
Ap, (= +u az ) = KyP( Ya Ya) + hy : (60) 
Conservation of moisture in the desiccant: 
1 fà] W hP Ta ai T, 
U E @ ieee gy es (61) 
Pia hii 


T.S. Ge et al. / Renewable and Sustainable Energy Reviews 12 (2008) 1485-1528 1503 


Conservation of energy in the process air: 


OTa OTa 
Apg(Cpa + Y aCpv) (T +u az 


) = hP(Ta = Ta) F KyPcpo( Ya = Ya)(Ta T Ta). 


(62) 


Conservation of energy in the desiccant: 


1 OTa 
> V alcpa + Wepi) +f sp€pspl Or 


= AP(T 4 = Ta) F K,P( Ya = Yast 
hP(T, — Ta) 
Nyy 


F KyPepy( Ya = YaXTa Ta) F dst: (63) 


This model was carried out to discuss the performance of RD silica gel desiccant wheel 
under high pressure. It has been reported that when the separation factor R is 0.01, 
operating pressure almost had no effect on desorption capacity, but sorption capacity 
enhanced with increasing pressure. Especially if the pressure is higher than 3atm, the 
sorption capacity could be greatly improved. When R reaches 1, the adsorption capacity is 
enhanced and desorption capacity gets worse with pressure rising. If the pressure is higher 
than 3atm, the variation of adsorption capacity is not significant. It was also found that 
desiccant wheel achieves a better performance at the condition of R = 1 as compared to 
R= 0.01. 

Condensation and frosting processes occurring in energy wheels have been investigated 
too. Simonson and Besant [28] expanded the model proposed by them [7] to include 
condensation and frosting processes. The difference between the governing equations of 
these two models is that the adsorption heat qs: in Eqs. (42) and (43) is replaced by heat of 
vaporization hy. In saturation condition, the density of water vapor can be approximated 
by 


Ps 


Sprene 4 
Po = Pos RT, (64) 


where R, is the specific gas constant of water vapor. 

The numerical schemes are similar to those used in Ref. [7]. Mathematically and 
physically the problem becomes more complex when condensation and frosting are 
included. Hence these conditions require a different numerical solution strategy and then 
different algorithms are adopted during sorption and saturation process as shown in 
Table 1. This model could simulate the condensation and frosting happening in energy 
wheel. The solutions were validated by comparing the simulated results with experimental 
data on a commercial energy wheel. This study analyzed the sensitivity of condensation 
and frosting to wheel speed and desiccant type. The simulation results showed that a 
desiccant with a linear sorption curve has smaller amounts of condensation/frosting 
compared with Type I sorption curve. 

Dimensionless equations of air-to-air energy wheel were derived using this model [29]. It 
was indicated that the performance of energy wheel is a function of the operating 
temperature and humidity. In addition, based on this model, Jeong and Stanley [30] 
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Table 1 
Different algorithms used to solve the governing equations [28] 


Sorption (®< 100%) Common Saturation (® = 100%) 


1. Estimate the rate of phase change and 
the properties needed in the governing 


equations 
2. Solve the Tm, Ps, W, u and T, 2. Solve the T,,,, m, W, u and T} 
fields in order with Eqs. fields in order with Eqs. 
(42)-(46) (42)-(46) 

3. Update the properties 
4. Update the rate of phase 4. Update the rate of phase 
change (m) change (p,) 


5. Return to step 2 and iterate until a 
converged solution is reached 

6. Adjust, as necessary, the algorithm 
used in each grid for the next time-step 
7. Increment time, return to step 2 and 
iterate until a quasi-steady solution is 
obtained 


Table 2 
Features of gas side resistance (GSR) model reviewed in present work 


Ref Dimension Condensation Boundary conditions Parameters of initial condisions Other relations 


Mode Expressions 


13] I = = - Eqs. (16), (25) 

14] I I Eqs. (5), (6) T, Y, Eqs. (26)-(30) 

6] I I Eqs. (5), (6) Ta Ya Ta Ya W Eqs. (26), (35)-(41) 
7 I Eqs. (5), (7), (55) - Eqs. (47)-(54) 

26] I a I Eqs. (5), (6) Ty Y, Eqs. (26)-(30) 

27] I af I Eqs. (5), (6) Ta Ya Ta Ya W Eqs. (26), (35)-(41) 
28] I J = Eqs. (5), (7), (55) - Eqs. (47)-(54), (64) 
8] n I Eqs. (5), (6) Ta Ya T4 W = 

Il Eq. (9) 


developed practical enthalpy wheel effectiveness correlations which can be more easily 
used to design and analyze passive wheel. 


3.2. Two-dimensional model 


Multi-dimensional mathematical models reflecting the effects in radial or circumferential 
directions could be used to analyze heat and mass transfer processes comprehensively with 
improved accuracy. But the complexity increases simultaneously. Usually two-dimensional 
model is recognized as a compromise method between accuracy and complexity. 

Yu et al. [8] have derived a two-dimensional GSR model. Convection items 0Y/d@, 
OW /Og and OT /ðọ in circumferential direction which reflect the heat and mass transfer 
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due to the rotation of the wheel were added into the model. The control volume of the 
model was the same as that shown in Fig. 2. The governing equations were as follows: 
Conservation of moisture in the air: 


OY, OY, m, OY, KyAy 
+o—*+ = 
Ot õp p;Fs Oz pF 


l 


(Ya Ya). (65) 


Conservation of moisture in the desiccant: 


ow OW _ Ky 
a Op pð 


( Ya Yq). (66) 


Conservation of energy in the air equation: 


OTa OTa m, OT, hAy 
D = 
Ot õp piFs z PiF (Cpa =P YaCpv) 


(Ta — Ta). (67) 


Conservation of energy in the desiccant: 


oT oT h K,(Va-— Y 
d d (T, — Ty) + A Das 
Ot 0p  paôlCpa + Wcpi) Pg(Cpa + Wop) 


(68) 


Boundary conditions of modes I and II were introduced to solve the problem. The 
governing equations were discretized into finite difference equation group. For space 
varying items, one order upwind difference scheme was adopted for convection terms. For 
time variation items, different schemes were used according to the demands of 
computation. To compute steady state, fully implicit scheme was employed because of 
its reliability of convergence. To investigate the dynamic behavior of the dehumidification 
device, explicit scheme which does not need the iteration process was recommended. 
Owing to the constraints of stability, time step should be small enough, or the divergence 
of solution cannot be avoided. Moreover, small time step was also numerically necessary 
for accurate computation. Crank—Nicolson scheme could be carried out for dynamic 
investigation too. A proper mesh size should be specified to save computation time. It was 
found that the predicted result of this model is more close to the actual process than that of 
Charoensupaya and Worek [13]. This is mainly because although the convection items in 
circumferential direction can be neglected under low rotation speed, they influence heat 
and mass transfer process in almost all occasions. Hence, when they are added to the 
model, more practical results can be achieved compared with one-dimensional model. Also 
the analytic results indicated that when the density of the adsorbent py is smaller than 0.2, 
the influences of OY /0g and OT /0@ cannot be negligible. 

Feng et al. [31] applied analogy method to solve the proposed governing equations by 
Yuet al. [8]. Their study have shown that the performance of dehumidifier is a function of 
the nondimensional structure parameters and thermophysical property parameters (such 
as nondimensional rotation period #. = (m,*t,)/(o,*Fs*L), nondimensional moisture 
ratio m = (p,° Fs* Ya)/(Wmax* My)) with the operation parameters keeping constant. The 
GSR models summarized in this part are reviewed in Table 2. 


4. Gas and solid-side resistance (GSSR) model 


The previously mentioned GSR models do not consider the heat conduction and mass 
diffusion within solid side, thus they do not reflect the actual transfer process occurring in 


Table 3 


Features of gas and solid side resistance (GSSR) models reviewed in present work 


Refs. Dimension Conduction in Diffusion in Boundary Parameters of Other relations 
solid-side solid-side conditions initial conditions 
Zz r @ 2 r Mode Expressions 
[32 I PGS model I Eqs. (5), (6) - Eqs. (26), (75), (76) 
II Eq. (9) 
[9] I PGS model - - - Eqs. (28), (83}-(88) 
[33 I PGS model II Eq. (10) Y, H, Eqs. (93)-(95) 
(«NTU Model) 
[34 I PGS model I Eqs. (5), (6) Ta Ya Eqs. (35), (96), (98), (105), (106) 
(Radial-flow wheel) 
II Eqs. (103), (104) 
[35 I J y — — — Eqs. (35), (96), (109), (112)-(114) 
[10 Il Ni J Ul Eqs. (11)}(14) Ta Ya Eqs. (119)-(121) 
[36 II J af af I Eqs. (5), (6) Ta Ya Ta W Eqs. (26), (36), (126)-(128) 
II Eq. (9) 
[37 II V J J af I Eqs. (5)-(8) — Eqs. (129), (132), (135)-(137) 
lll Eqs. (12), (14) 
[38 1/11/01 J af aff J I Eqs. (7), (8) - Eqs. (140)-(142) 
lll Eqs. (12), (143)-(145) 
[39 I PCP model I Eqs. (153), (156) W T, Eqs. (36), (126)-(128), (158)-(161) 
(Diffusion in the radial mM Eqs. (154) 
direction of a particle) 
[40 I PCP model H Eqs. (166), (167) W H, Ha Ya Eqs. (162), (164), (165) 


90ST 


8ZSI-S8FT (8007) ZI SMamay Absauq ajqouipjsng pun ajqomauay / 7D Ja aD `Ç iL 
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the desiccant wheel. To solve this problem, models considering heat and mass transfer in 
both solid and air side are proposed. Such models are called GSSR models (listed in 
Table 3). GSSR model is more preferable compared with GSR model for it can explain the 
effect of heat and mass diffusion within solid desiccant and the precision of governing 
equations describing the physics of desiccant wheel is enhanced. 

Typical diffusion theory holds that water vapor molecule has three types of diffusion 
mode in porous mediums which includes: ordinary diffusion, Knudson diffusion and 
surface diffusion. The actual diffusion process is the combination of these three 
mechanisms. 

Ordinary diffusion which is also called molecular diffusion or Fickian diffusion 
occurs when the mean free path is relatively short compared to the pore size. This 
diffusion mode is applicable to Brownian motion, where the movement of each 
particle is random and not dependent on its previous motion. The transfer diffusivity 
Do relates the macroscopic flux of molecules in a system to a driving force in the 
concentration. Knudson diffusion occurs when the mean free path of gas molecules 
is relatively long compared to the pore size. It describes situations in which gas mole- 
cules collide more frequently with flow boundaries than with other gas molecules. 
Knudsen diffusion is dominant for pores that range in diameter between 2 and 50 nm. 
Surface diffusion is also used to explain a type of pore diffusion in which solutes 
adsorb on the surface of the pore and hop from one site to another through inte- 
ractions between the surface and molecules. All of these transport processes can be 
described by means of Fick’s law of diffusion with appropriate diffusion coefficients 
(Do, Dx, Ds). 

Pesaran and Mills [39] investigated the diffusion mechanism in silica gel. It has been 
reported that surface diffusion is the dominant one for moisture transfer in micro porous 
silica gel such as RD gel (regular density silica gel). While for macro porous silica gel, for 
example ID gel (intermediate density silica gel), both Knudsen and surface diffusions are 
important mechanism. Ordinary diffusion can be negligible for the silica gel at atmospheric 
pressure. 

Based on different methods to account for the SSR, GSSR models are classified into 
three subcategories: (i) pseudo-gas-side (PGS) models in which lumped convective 
coefficients are adopted to account for SSR; (ii) gas and solid-side (GSS) models in 
which the second order diffusion terms are added to account for SSR; (iii) parabolic 
concentration profile (PCP) models which assume that a parabolic concentration profile 
for moist concentration exists at all the times in the desiccant particle. 


4.1. Pseudo-gas-side (PGS) model 


A simplified method to include the heat and mass transfer in the solid desiccant is using 
lumped transfer coefficients (h' as heat transfer coefficient and K, as mass transfer 
coefficient) which include heat and mass transfer resistances not only in gas side but also in 
solid side. Such modified models are named as PGS model. Since the second-order 
diffusion terms are not added, the complexity and computational time of PGS model do 
not increase compared with GSR model. The lumped transfer coefficients are usually 
obtained by the correlations of experimental data. The correlations for silica gel 
as adsorbent given by the Hougen and Marshall [41] based on Ahlberg’s experimental 
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data [42] are in wide use: 


h = 0.683m, Re~” cpa, (69) 


K, = 0.704m; Re?! (70) 


4.1.1. Axial-flow PGS model 

Maclaine-Cross [43] has presented a finite-difference model known as MOSHMX 
(method of solving heat and mass exchangers). The gas-side controlled, lumped 
capacitance mass transfer coefficient was used in the model which is actually a typical 
PGS model. A model of desiccant wheel with PGS coefficients similar to the one proposed 
by Maclaine-Cross [43] has been described by Holmberg [32]. In the model, air thermal and 
moisture storage capacities OY /Ot and OT /Ot were neglected due to their small value 
compared to the matrix capacities. The control volume is depicted in Fig. 7. 

Mass and energy conservation equations were: 


OYa , fittmOW _ 


Meee u oy 

rita oe J (= tin kAm a) =0. (72) 
Mass and energy transfer equations were: 

ee eV), (73) 

ie oH _ a | x (Ta = T) F hyl n= Fal (14) 


where it is assumed that the rotor matrix is composed of supporters (sp) and adsorbent (d) with 
the mass fractions 1—f and f, respectively. Thus the enthalpy of the matrix is defined by 


Hm = (l —f\H sy +flHa TF W(H ws — fst T hyy)). (75) 
The equilibrium isotherm was in the form of Eq. (18) 
W = 0.558. (76) 


Process air 


. Poe OSS Sa 7 
rotation | Aie 
l — 
Regeneration air i Tas Yg 
= 
l 
yA | Heat | Mass 

se 
T 
l Tm W Z 
i Matrix o 


Fig. 7. The model of control volume and a side view of one of ducts for axial-flow PGS model [32]. 
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Holmberg adopted implicit Crank—Nicholson scheme for the energy conservation 
equation (72) and ordinary central difference schemes for other equations. Gauss-Seidel 
solution technique was employed to solve the model. A parametric study was conducted to 
show the influences of matrix heat and moisture capacities on temperature and humidity 
efficiencies. It was indicated that there is an optimum drying humidity efficiency with 
variable heat capacities and the optimum value increases in magnitude with increasing 
moisture capacity. 

Based on the studies of Maclaine-Cross [43] and Holmberg [32], Jurinak and Mitchell [9] 
neglected the axial heat conduction inside the matrix and developed a dimensionless PGS 
model in period j (j is a period index) of rotary dehumidifier. The control volume was the 
same as shown in Fig. 7 and the governing equations were given as: 

Mass and energy conservation in period j: 


OYa. OW _ 
0H, y OH in _ 
aor t GU Zar = 0: (78) 


where coefficients u; and x; are defined as: 
H; = Mary matrix mass, j/ Madry air containned in matrix, j kj = L/(u;9)). (79) 


Mass and heat transfer equations in period j: 


Ya 

az = Awj(Ym T Ya), (80) 

OH, 

a = Ay j(LeoCpg(Tm a Ta) F H,( Yin = Ya)), (81) 
where dimensionless coordinates in the governing equations are 

z =z/L, 0 =0/t, 0; = 0;/t, (82) 


and dimensionless heat and mass transfer length are given by 
4h' 


Arj =Jrj;L[v, GEF (83) 
Awj = JwjL/v, Jwj = ‘Ky : (84) 
Pade 
Overall Lewis number was expressed as 
Leo = Jrj/Jwj. (85) 


The enthalpy of the matrix at any temperature and water content with respect to a dry 
matrix at 273 K were given by 


wW 
Hm = Com(Tm — 273) + cC W (Tm — 273) + hg | (=) dW, (86) 
0 


it pa eye eS 


A= tant 2 
ia =en : constant, (28) 
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where the dry matrix thermal capacities, Cpm, is related to the specific heat of desiccant cpa 
and of the nonsorbing matrix materials c,; 


Com = Cpd F DLomepd/f- (87) 


Six different isotherm equations were adopted in the study to investigate their influences. 
They are functions of W/W max and here only the simplest linear style is listed in detail: 


P= G(w/ W max) = W/ W max- (88) 


A finite difference solution was used to study the effects of different matrix properties on 
the performance of a rotary dehumidifier, and in each case the results were obtained for the 
maximum dehumidification efficiency. The results were useful to select desiccant and 
analyze the open-cycle systems. 

Cejudo et al. [44] validated the model proposed by Jurinak and Mitchell [9] by 
experimental data. Moreover, they used the real data to train a neural network to calculate 
the parameters of a new neural network model. Their results demonstrated that the new 
neural network model is more adequate to simulate wheel performance when heat losses of 
the rotating wheel are taken into account. 

Later, Van den Bulck et al. [33] pointed out that, the model solution in Ref. [9] involved 
considerable computation and was not suitable for long-term simulation of desiccant 
cooling systems. In order to reduce the calculation efforts, they introduced number of 
transfer units (NTU) in the governing equations. The new conservation equations using 
NTU for period j were 


a oF = ce 
a ie “= 0, (90) 
“ = NTUy,(Ymn — Ya), (91) 
— = NTU, rar n — Ta) + NTU yj Hol Ym — Ya), (92) 


where air stream mass flow rates, wheel rotation speed and mass of desiccant are combined 
into two dimensionless capacitance rate parameters T;: 
Med 


f=— 
ttn, 


0 1 1 
i= = —— Llen amm 
G=1,2), t T (o<r<7), Z=z/L (93) 


and the number of transfer units is defined as 


NTU 2 94 

wj — ; > ( ) 
mj 
h A; OH, 

NTUy =o u= (See) e 
ljCpg g/ ¥, 
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This approach is called effectiveness-number of transfer units (e-NTU) method which 
clearly represents the influence of flow unbalance and transfer coefficients respectively. 
Boundary conditions in mode II and thermophysical property parameters obtained from 
Ref. [45] were applied to close the partial difference equations. Then the model was solved 
by means of characteristics and the shock wave method. Using this model, the 
performance correlation of ideal dehumidifier with infinite heat and mass transfer 
coefficients was calculated as a function of NTU and operating parameters T. 

Subsequently, correlations for the humidity and enthalpy effectiveness of actual silica gel 
dehumidifiers with finite transfer coefficients as functions of NTU were conducted by Van 
den Bulck et al. [46]. In their analysis, equilibrium isotherm and enthalpies of gas and 
matrix were calculated by 


P = (2.112 W)*/" (29.91 Py)" (96) 
Hy = Cpg(Tg — Tret) + Ail Tret] Yg, (97) 
WwW 
AEN = Tae J lT iet] — Gul, Tret) 42. (98) 
0 


The correlations were validated by comparing the effectiveness values obtained from 
these correlations with effectiveness values provided by MOSHMX [43] and the ideal 
dehumidifier theory [33]. The results showed that a e-NTU model incorporating these 
correlations can rapidly and accurately predict the dehumidifier performance. And they 
were suitable to be used in long term performance simulation of desiccant cooling system. 

Also, by using the same e-NTU model, Klein et al. [1] analyzed a regenerative enthalpy 
exchanger, Van den Bulck et al. [47] investigated a desiccant cooling system and a modified 
method was proposed to solve s-NTU model for the desiccant which has discontinuous 
isotherms and isopiestics such as LiCl [48]. 


4.1.2. Radial-flow PGS model 

The one-dimensional desiccant wheels focused on axial flow have been extensively 
investigated in the literatures as described in previous section. Elsayed and Chamkha [34] 
studied the radial flow desiccant wheel and developed a model to predict its performance. 
The radial flow wheel with an inner radius rı and outer radius r2 is shown in Fig. 8. It is 


Regeneration Dehumidified 


Fig. 8. Schematic of radial-flow desiccant wheel and the model of control volume for radial-flow PGS model [34]. 
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divided into two sections: period 1 of fraction f for regeneration air and period 2 of 
fraction 1—f for process air. The air flows in the radial direction inwards or outwards. The 
assumptions for the model were similar to that of Jurinak and Mitchell [9] and Van den 
Bulck et al. [46]. 

Conservation of moisture in the air was expressed by 


Aj0Ya ; 
res dr = K’ Ay dA; dr(Y¥m — Ya). (99) 
A; or J 
Conservation of moisture in the matrix was given by 
2nrdr OW 
= dt = K',Ay dA; dr(Ya— Ym). (100) 


Mn®—>— —— 
"mrs — r?) Ot 


Conservation of energy in the air yielded 


d4; 0H, j 
ti ar + Ky Ay dA; rH (Ym — Ya) = hAy dA; dr(Tm — Ta). (101) 
i or > 
Conservation of energy in the matrix yielded 
2nrdr Hm . dA; 0H; 
Mmo — L dt + mi 1 dr =Q. (102) 
n(r3 —r{) ôt A; Or 


The boundary conditions in modes I and II were applied in the model. Periodic 
boundary condition was expressed by the following equations: 


W(r/ri ot = 0) = W(r/r1, ot = 1), (103) 


Tm(r/ri ot = 0) = T,,(r/ri, ot = 1). (104) 


Based on experimental data, the corrected compound heat and mass transfer coefficients h’ 
and K’ were dependent on the superficial air velocity u within a matrix cross-sectional 
area [49] 


W = Cy’, (105) 


Ki, = Cmt pas (106) 


where Ch, Cm and « are constants. 

The relations to determine the heat of sorption, equilibrium humidity ratio of air in 
contact with a silica gel matrix and the enthalpy of the silica gel were the same as Eqs. (35), 
(96)-(98). 

These differential governing equations were reduced to dimensionless NTU style and 
discretized by finite difference technique on staggered grid. Parametric studies were 
performed using the model to determine the effects of flow configuration and operational 
parameters on the performance of radial flow silica gel desiccant wheel were studied. Also 
axial flow desiccant wheel was compared with conventional radial flow wheel. It was 
founded that whose performance is better depending on the values of the design and 
operating parameters. For instance, in ventilation and recirculation modes, a radial flow 
wheel with « = 2 (Eq. (105)) has lower outlet humidity ratio of process air and higher 
performance coefficient than axial flow wheel. But when «>6, the results are opposite. 

It should be noted that composite transfer coefficients h’ and K’, are used to account for 
solid diffusion resistance in PGS models. Therefore, although second-order differential 
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terms are not added in the governing equations and the style of the equation is almost the 
same as GSR model. The SSRs are indirectly included. That is why PGS model is excluded 
from GSR model and included in GSSR model. 

Since the composite coefficients are obtained from large numbers of experimental 
results, it is not convenient for studying a variety of desiccant materials and geometries. 
And by now, only the data for silica gel has been obtained. Hence, the reliance on such 
excessive available data limits the application of PGS model. Besides, the reliability of PGS 
model is not good enough which hinders its further development. Clark et al. [50] tested a 
prototype bed designed for solar air conditioning and found the results were in rather poor 
agreement with predictions based on the PGS model, particularly after a step change in 
inlet air condition. They concluded that the discrepancy was due to the shortcomings of the 
model, since the SSR exceeds the GSR under these conditions. The experimental result 
reported by Clark et al. [50] was somewhat limited and imprecise for the bed with a large 
prototype design. Hence Pesaran [21] performed bench scale experiments on thin beds of 
regular density silica gel with step changes in inlet air humidity. Their results conformed 
that the SSR was indeed generally larger than GSR. Also, the results of Mei and Lavan 
[51] indicated that the predicted results obtained from PGS model were acceptable 
compared with experiment results, but the coherence was not very well. Furthermore, 
Gheelayagh and Gidaspow [52] pointed out that the composite coefficients h’ and K’, that 
represent the SSRs were not constant, their value changed with the time. So adopting 
constant composite coefficients cannot reliably forecast the performance of dehumidifier. 


4.2. Gas and solid-side (GSS) model 


Since PGS models adopting lumped transfer coefficients cannot provide very accurate 
results, and there exists a demand for suitable models containing SSRs. With the increasing 
knowledge of diffusion mechanism inside desiccant, gas and solid-side (GSS) models are 
proposed in which heat and mass transfer diffusional terms in solid-side are directly added. 
Although GSS model needs more time to compute the additional second order diffusion 
terms compared with PGS model, it is more related to the true physics of the problem. 
Therefore, GSS model is likely to give better agreement with experiment results. 


4.2.1. One-dimensional model 
Using the control volume illustrated in Fig. 9, San and Hsiau [35] have derived a one- 
dimensional model considering the axial mass diffusion and heat conduction in solid 


Fig. 9. The model of control volume for a one-dimensional GSS model [35]. 
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desiccant, while the storage terms of moisture and energy in air O0Y,/0t and 07 ,/0t were 
neglected. The mass and energy balances were expressed as: 
Conservation of moisture in air: 
OY, AyV 
Maaz OL 
Conservation of moisture in solid: 


Ow 5. ow 
Wal ae AVG ga 


where the surface diffusivity is evaluated by 


K,( Yn T Ya). (107) 


= Ay VK,( Ya = Ym), (108) 


1.6 x 10% 
Ds = 2” - Pm exp( 0.974 x 1078 x a), (109) 
Conservation of energy in air: 
. ôT, AyV 
Cpgg ap OL A(T m — Tq). (110) 
Conservation of energy in solid: 
ôT n ô Tn 
MmCpm ar AyV 5k a2 = Ay Vh(Ta — Tm) + Ay VK,y( Ya — Ym)Qst (111) 


where k is the compound conduction coefficient and calculated using the following 
expression: 


k=f-ka + —fkep (112) 


taat 


and subscript “m” means the solid part of the wheel which contains desiccant and 
substrate materials and then “m,” stands for the total mass of solid desiccant wheel. 

RD gel was selected as desiccant in the study. Due to surface diffusion is the only 
mechanism need to be considered for RD gel [39], only “Ds” which represents surface mass 
diffusion coefficient appears in the model. The heat of sorption and equilibrium isotherm 
relationship were determined by Eqs. (35) and (96). The saturation pressure of vapor P,, 
was evaluated using the following expressions: 


Ps z a+bz+ cz? 
log10 318.1677 T IFd ia 
z = 100.97 — T; a= 3.2437814; b= 5.86826 x 1073, 
c = 1.1702379 x 1078; d = 2.1878463 x 107°. 
The relationship between Y,,, Tm and ® was obtained from ASHRAE data 
3 
Y,, = 6.22 x 10°@ (114) 


10421429-17.5(Tm—273.15)/(Tm—=35.85)] = ®/100 i 


A mixed numerical scheme was developed to solve the governing equations. The two 
governing equations in desiccant side (Eqs. (108) and (111)) were solved using full-implicit 
finite difference scheme and the other two in gas side (Eqs. (107) and (110)) were solved 
using backward finite difference scheme. The full-implicit finite-difference scheme was 
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verified to be effective for solving the equation with diffusion terms. Also this mixed 
scheme was compared with the explicit scheme which solved the four governing equations 
totally with a backward finite difference scheme. It showed that for the same accuracy the 
mixed scheme takes less computation time than the explicit scheme. On the other hand, for 
the same time step At the balance in the mixed scheme is much better than that in the 
explicit scheme. 

Numerical results of this model indicated that the axial heat and mass diffusion are 
governed by two important parameters NTU and 67/4Bi, and for a case of zero SSR, the 
calculated performance of the silica-gel desiccant wheel is poor. Furthermore, it was found 
that only taking account of surface diffusion is good enough for predicting RD silica gels 
desiccant wheel, but for ID silica gels the Knudsen and ordinary diffusion are needed to be 
considered. 


4.2.2. Two-dimensional model 

Charoensupaya and Worek [10] established a two-dimensional GSS model in which the 
heat conduction and mass diffusion including gas phase diffusion and surface diffusion in 
the radial direction within solid desiccant were considered. The control volume is 
illustrated in Fig. 10. 

Conservation of moisture in the process air was expressed as 


ha [10Ya  O¥g 
7 ( + ) = 2K (a= Ye) (115) 


u Ot Oz 
Conservation of moisture in the desiccant felt was given by 


dYa ow Ya ow 
€tPad Ot +( E1)Pa Or DGPad Ore Dspq Or2 = 0, (116) 


where Dg and Ds represent the effective gas phase diffusivity and surface diffusivity, 
respectively. 
Conservation of energy within the desiccant felt gave 


OH ad aTa Ta Ow 
& Pad a; + (1 — &)paCpa ar ka a = q! E)Pals ar> (117) 


where the subscript “‘ad’’ means the air in desiccant pore. 


nonhyroscopic wall 


Yv 
rA m w 


Sorbent t ag 
Adsorption air — ______ Desorption air 
=> m ae aa abs aise aq <i Sul 
Ta Ya 
S“ v 
A t 
ee ee 4 L. z 


Fig. 10. The model of control volume for a two-dimensional GSS model [10]. 
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The rate of energy transfer between process air and desiccant felt yielded 


ôH, ðH, 2K, 
+u = 
Ot Oz Pga 


OH, 2h 
(Yalt, x, %a) — Ya) a + (Talt, x,%a) — Ta). (118) 


gra 


The boundary conditions were in the form of mode II. Also, the temperature of the wall 
was assumed to be the same as the desiccant temperature at r = 0 


T(t, z) = Ta(t, z, 0). (119) 
The model was verified using the following equilibrium relation and heat of sorption: 

W = 0.59, (120) 

Ist = hy. (121) 


The heat and mass transfer coefficients were approximately by constant value that can be 
calculated from the Nusselt and Sherwood number correlations. 

An explicit finite-difference scheme was adopted as the solution method to calculate the 
model. The model was used to analyze the effect of some parameters, including heat and 
mass transfer Biot number, on the performance of desiccant cooling system. It was 
revealed that there exist optimum values of COP and cooling capacity for given operating 
conditions. The optimum values can be reached by controlling the system cycle time. Also 
the accuracy of this model has been verified by experimental data in Ref. [53]. 

Yu et al. [36,54] also constructed a two-dimensional mathematical model for desiccant 
wheel based on the control volume shown in Fig. 2. Both the mass diffusion and heat 
conduction inside the solid desiccant wheel in axial as well as circumferential directions 
were considered in the model. The governing equations are listed as follows: 

Conservation of moisture in the air: 


OY, OY, m, OY, K,Ay 
D + = 


= Ya— Yq). 122 
a tO Gp + pF, a pf TT (122) 
Conservation of moisture in solid adsorbent: 
ow oW 2ln(r2/r)ð® W CW) K,Ay 
DU — F, =— Y,— Yq). 123 
Pap | oe apt? ae |, en (123) 
Conservation of energy in the air: 
OTa OTa m, OT, hAy 
: = Ta — Ta). 124 
ot iki dp T piFs 0z pF (Cpa + Yacn) ) oe 


Conservation of energy in solid desiccant: 


OTa P oTa kal — F;) 2ln(r2/r1)ð Ta OT 4 
Or Op — [PalCpa + Wepi) + PspCpsp] | r3- r? 8G? ~— 82? 
— [hA v(Ta = Ta) + K,A v( Yq = Yaastl 
E [Pa(Cpa F Wep1) T Psp€psp] l 


(125) 


The boundary and initial conditions were similar to those of Ref. [14]. Eq. (26) was 
adopted to calculate Y4. Equilibrium isotherms ®, obtained by fitting fourth-degree 
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polynomials to the manufacture data [21], and heat of sorption qs, for RD and ID were: 
For RD gel: 


® = 0.0078 — 0.05759 W + 24.16554W> — 124.478W? + 204.226", (36) 
—12400W +3500 W<0.05 (kJ/kg water), 
ds = $ _1400W +2950 W>0.05 (kJ/kg water), (120) 
For ID gel: 
_ f 1.235W + 267.99W° — 3170.7W° + 10087.16W* W<0.07, 127) 
~ ) 0.3316 + 3.18 W W>0.07, 
—300W + 2095 W<0.15 (kJ/kg water), 
Iss = $ 2050 W>0.15 (kJ/kg water). 28) 


Central difference scheme was applied to discretize the diffusion items and the same 
schemes as in Ref. [8] were adopted to discretize the other terms. It was demonstrated the 
numerical results of this model are very close to the steady and transient experimental data 
of Refs. [55,56]. 

Dai et al. [57] pointed out that the model proposed by Yu et al. [36] was helpful for 
predicting the performance of desiccant wheel but was not brief enough. Hence they 
adopted psychrometric wave chart analysis associated with the finite difference model 
trying to obtain the numerical results efficiently. In their study, RD silica gel was used as 
adsorbent. According to the conclusion of Pesaran [39], the surface diffusion is dominant 
for RD silica gel, and then the effective diffusivity was given as 


De = Ds = Do exp (-o.074 x 1076 x t), Do = 0.8 x 107°. (129) 
d 

Using the wave analysis in psychrometric chart incorporating with numerical 
calculation, the influences of some important parameters, such as heat capacity, 
adsorption heat, rotation speed, thickness of the desiccant matrix, etc., on the performance 
of RD silica gel desiccant wheel system were discussed. It was found that the chart method 
is rapid in evaluating the system performance and is useful in system design. 

Gao et al. [58] have described a mathematical model that can be seen as the 
simplification form proposed by Yu [36], where diffusion terms in circumferential direction 
were ignored and the sensible heat transfer was included in the governing equations. They 
conducted a parameter analysis employing the model to investigate the effects of passage 
shape and geometrical size of the matrix on the performance of desiccant wheel. 

Zhang and Niu [37] have presented a two-dimensional transient model in which heat 
conduction, surface and gaseous diffusion in axial as well as radial directions of the 
desiccant wheel are included simultaneously. Due to symmetry, the mid-plane of a channel 
was treated as adiabatic wall, and a half-size channel was chosen as the control volume as 
shown in Fig. 11. The derived governing equations were: 

Conservation of moisture in the air stream: 


10Y, OY, 4K, 
= Y,— Ya). 1 
u Ot s 0z Tau ) (a 


1518 T.S. Ge et al. / Renewable and Sustainable Energy Reviews 12 (2008) 1485-1528 


C Air Stream 


Fig. 11. The model of control volume for a two-dimensional GSS model [37]. 


Conservation of moisture in the desiccant: 


Ya Ow ð Ya 0 Ya 
Pad Sy + Pd Bp = Pad = (2 Oz ) Te, (2 or )| 


0 Ow 0 Ow 
+ pals (0s z) ty (0s or J|; a) 


where the effective surface diffusivity Ds can be calculated by Eq. (129) and Dg which 
stands for the composite effective diffusivity of Kundsen and ordinary diffusion is given by 


e/a \ Pee ue iy 
Dg = — | — + — Do = 1.75 10° “—, Dx = 97a| — 132 
G (Et) š o 758 x 10 B’ K Ta m 5 (132) 
where a (m) is the pore radius of the adsorbent and mı (kg/mol) is the molecule weight of 
water. 
Conservation of energy in the air stream: 
10T, OT, 4h 


je ae (133) 


Conservation of energy in the desiccant: 


aTa Ta Ta Ow 
Pit a, = i( ga t apa) TIP ae (134) 


where c,, is the total heat capacity of wet desiccant, which included two parts: dry 
desiccant and adsorbed water and is calculated by 


Cpt = Cpa + Wepy. (135) 

By using heat mass transfer analogy, the relations between h and K, can be expressed as 
h 

K, = —. 136 

? CpgLe (136) 


Water content in the desiccant is governed by a general sorption isotherm as 
(137) 


The boundary conditions were similar to Eqs. (5)-(8), (12), (14). This model was 
normalized and numerically solved by means of alternating direction implicit (ADI) 
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method [59]. Using the model, effects of rotary speed, NTU and specific area on the 
performance of silica gel wheel were studied and compared in desiccant wheels for 
dehumidification and enthalpy recovery. 

Niu and Zhang [60] performed an analysis study adopting the above model to study the 
influence of wall thickness on the optimum rotary speeds for dehumidification desiccant 
wheel and enthalpy recovery wheel. 

Sphaier and Worek [38] derived a desiccant wheel mathematical model which 
accounts for heat and mass transfer processes inside solid desiccant too. The model is 
different from others in two aspects: (1) the diffusion terms in the governing equations 
are given in the form of vector, hence governing equations can be easily converted 
into one or two, even three-dimensional ones; (2) parameters describing detailed 
constitution of felt are included in the model. The model is thought to be useful for 
analyzing the transport phenomena within the porous adsorbent felt and the influence of 
supporting structure. It is assumed that the entire felt (subscript “‘f’) is composed of rigid 
matrix (subscript ‘‘fm’’) and pores (subscript “‘fp”) in which both the gas and adsorbed 
liquid phase coexist. The matrix includes solid adsorbent desiccant (subscript ‘‘fd’’) and 
inert materials which has no effects on the sorption characteristics of adsorbent and 
surface diffusion. 

Based on the control volume shown in Fig. 12, the derived conservation equations of 
moisture and energy in the air stream were similar to Eqs. (130), (133), and conservation 
equations of moisture and energy in the desiccant were: 


Yp OW fim 

rpa ae + Eppa a = paV (Deen Y Yp) + pef VDs er VW tin) (138) 
oT © W m 

prepa = V{ky er VT) + paf | — £) ar — V(Ds ett V W jin) | dst» (139) 


where the effective heat conductivity kyr gas phase and surface diffusivity Dg orp, Ds cer are 
defined as 


kf ett = Etkfp + (1 — £1)kfn, (140) 
DGett = &DG/Tg, (141) 
Ds et = (1 — &)Ds/Ts. (142) 


Eqs. (7) and (8) were utilized to depict the insulated and impermeable boundaries. 
Boundary conditions in mode III were obtained through interface balances: at r= Rg 
(Fig. 12), the boundary condition for mass transport was the same as Eq. (12) and for heat 


Impermeable wall, w 


A r =R, 
Sorbent felt, f r=Ry 
T r=R, 
| 
Process stream, g he i 
y Ji < | 


L. 


Fig. 12. The model of control volume for a GSS model with vector form [38]. 
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transport was given as: 


OT 0 m 

—ky eff e = NT _ Ts) F Ky(Va Eg Yp)dar T Pfa Ds,eft dst (143) 
and at the wall the boundary conditions were 

—PaD Gett Se = Pfa Ds,eff Or = 0, (144) 

: oT oY OW fin 

qs = -kjen A — p DGert 2 dar — Pya Ds,eff = Hi, (145) 
where qar is termed the heat of sorbate transfer at the interface 

4AT = Cp Ta = Ta). (146) 


Also the selecting equilibrium isotherm was in the form of Eq. (17). 

These governing equations were fully normalized using classical dimensionless groups 
for heat and mass transfer, and a numerical algorithm based upon the finite volumes 
method was employed to solve it. It showed that simulation results present a very 
reasonable agreement with the experimental data. Moreover the simulation results of a 
test-case were presented, suggesting a possible optimization to wheel initial-cost and 
compactness by reducing felt thickness. This model can be applied for the simulation of 
several configurations such as desiccant wheels or enthalpy exchangers. 


4.3. Parabolic concentration profile (PCP) model 


Pseudo-gas-side (PGS) model and gas and solid-side (GSS) model have their own 
shortcomings. The GSS model requires greater computational effort than PGS model 
because of the additional second-order heat and mass transfer diffusional items. However, 
the PGS model requires extensive experimental data for determination of the lumped 
transfer coefficients with different desiccant materials and its reliability is not good 
enough. Hence, a computationally efficient model that effectively considers both the SSRs 
and GSRs in desiccant wheel is needed. Parabolic concentration profile (PCP) model 
assuming an analytically convection concentration profile within the desiccant particle to 
account for solid-side diffusion is developed to achieve this objective. Usually, based on 
the assumed concentration profile, a differential governing equation describing mass 


12 


\ N H, T, Gas phase 


aa 


h(T;-T,) 


PE a ee, 


i 


I" 


Fig. 13. Schematic figure of: (a) Diffusion through a spherical particle. (b) A side view of one of ducts for a PCP 
model [39]. 
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diffusion within a desiccant particle is developed and included in PCP model to represent 
the moisture conservation in solid desiccant. 

Pesaran and Mills [39,61] proposed a PCP model for simultaneous heat and mass 
transfer in a thin packed bed of silica gel particles. It believed that diffusion of moisture 
through the porous particle consists of Knudsen diffusion and surface diffusion. The 
spherical particle model and adsorption bed model are illustrated in Fig. 13. 

The general form of moisture conservation equation for the spherical particle with radial 
symmetry was derived and given by: 


ow 10 > OW 
oe wee (2er F) In 


where D, is the total effective diffusivity which includes Knudsen diffusion and surface 
diffusion: 


g(W) 


P 


De = Dstt + Dk,eff 


oY 
í = m 14 
> g(W) (Sw). (148) 
For the thin beds, storage term 0Y,/0t was negligible and moisture conservation equation 
in the air was 
OY, 
g 5 £—K,(Y,— Y, — Y,)P. (149) 


z 


Conservation of energy in gas-side was described by 
. OT, 
Cogg 5> = —P[h + cw: K (Ys — Yg)]- (Tg — Ts), (150) 


where axial and radial conduction and the storage term 07',/0z have been neglected and 
the bed is assumed to be adiabatic. 
Energy conservation in the solid desiccant yielded: 


OT, 
ar = PMT, — T) - aK (Ps — Yh (151) 


where subscript “m” represents the adsorption bed. 
Eqs. (147) and (149) were linked by the equilibrium relation at the particle surface 


A Pm Com 


Y (z, =([W(r = R,z, t), T(z, 0), Ps]. (152) 
The boundary conditions and initial condition for Eq. (147) were 

mn =o, (153) 

pare r=R = Ky[Y5(z,0 — Y,(z, Ð], (154) 

W(z,r,t = 0) = Wo. (155) 


The boundary condition for Eq. (150) was 
T(z = 0,0) = Tin. (156) 
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While the initial condition for Eq. (151) was 


T;(z,t = 0) = To. (157) 
Heat and mass transfer coefficients were obtained by experimental correlations [21]: 

K, = 1.70m, Re”, (158) 

h = 1.60m, Re cg. (159) 


The properties for RD gel and ID gel were the same as Eqs. (36), (126)-(128), and specific 
heats of the gas and felt were 


Cpg = 1884Y, + 1004(1 — Y,), (160) 


R 
Com = 4186W avg + 921 Wavg = f 4nr W pp dr/(4/3ynR? pp. (161) 


The governing equations were put into dimensionless form and solved by finite 
difference method. Crank—Nicholson method was adopted for Eq. (147). Implicit Euler 
method was applied for Eq. (151). A fourth-order Runge-Kutta technique was used for 
Eqs. (149), (150). The transient results of this model were compared with the predictions of 
PGS model. It was proved that the simulated results getting from the new model agree 
better with the experimental data than those from PGS model for both ID gel and RD gel 
desiccant material. 

Chant and Jeter [40] performed a steady state and transient analysis of rotary desiccant 
dehumidifier using a PCP model. The heat transfer equation and conservation equation in 
the desiccant side were similar to those of PGS model (Eqs. (77), (78), (81)), except for gas- 
side transfer coefficients h and K, not compound transfer coefficients like h’ and K’, were 
adopted in PCP model. For mass transfer analysis, a parabolic concentration profile was 
assumed to exist in desiccant particle at all times, and for a spherical particle the profile 
was 


W =a, +a(r/R). (162) 
The moisture conservation equation for a spherical particle was 
ow 3K, 
ar = YW, Ta) — Yal, (163) 
Ot Pparticle” 


where W, is desiccant surface concentration which is derived using Eq. (162), Y4 stands for 
moist air equilibrium humidity ratio, they could be calculated by the following equation: 


KP parti 

W, = — Pe TY (Ws, Ta) — Yal + W, (164) 
SPparticlee 

Ya=ootqWt+oHgt+ Yq. (165) 


[62] The enthalpy and equilibrium adsorption isotherm expressions were similar to 
Eqs. (97)-(98) and Eq. (17). 

First the governing equations (Eqs. (77)-(78), (81), (163)) were associated with the inlet 
states of moist air for processing and regeneration as well as the desiccant bed initial state 
to construct transient model, which was used to study startup effects. Then periodic steady 
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state model was distinguished by applying the following periodic boundary conditions: 
Ha(z, t = 0) = H4(z, t = tr), (166) 


W(z,t = 0) = W(z,t = t,). (167) 


Transient model and periodic steady state model constituted the whole PCP model. The 
transient model was solved using Bulirsch-Stoer integrator method [63]. An implicit 
periodic steady-state solution with a linear adsorption isotherm (Eq. (165)) which utilizes a 
sparse matrix solving routine was developed to solve periodic steady-state model. 

Compared with the transient experimental results of a fixed bed dehumidifier. It was 
found that the PCP model does not predict well in the initial adsorption period, which is 
due to the reason that the moisture concentration profile is not parabolic until later in the 
sorption period. But the simulation results become similar to that of the PGS model as 
time moves on. Also the periodic steady state of a desiccant wheel system predicted by this 
numerical model has shown good agreement with the testing data of practical desiccant 
wheel at low and medium rotation rates condition. The authors discussed the efficiency 
and flexibility of PCP model too. Their results indicated that the model only requires 
1-3 min on a 25 MHz 80486DX PC for the cased considered. And since the widely used 
adsorbent property effective diffusivity De is used in mass transfer equation, it has the 
capability to simulate a variety of desiccant materials. Hence, PCP model is extremely 
efficient and an excellent tool for investigating alternative desiccant materials. 


5. Empirical model 


Beccali et al. [64] believed the models based on the detailed heat and mass transfer 
phenomena inside desiccant wheel are too complex. The solving process is complicated and 
time-consuming. Hence, in order to evaluate the performance of various kinds of rotary 
desiccant wheels available in the market, they analyzed the available measurements data 
from the manufacturer thoroughly and developed empirical model for evaluating the 
performance of silica gel desiccant wheel. It should be emphasized that the units of 
enthalpy, temperature and humidity ratio are KJ/kg, °C and g/kg, respectively, in the 
study. The derived correlations were expressed as 


Yout/Tout = (Fy Y? Ea Fy Yin + F3)Tin + (F4 Y, F Fs Yin F F6), (168) 
where F; is a second-order polynomial of Teg and can be expressed as follows: 


i = (A;T2,,+ BiT eg +C), i=1,...,6, (169) 


reg 


where 4; B;, C; are second-order polynomials of X,eg» 


Ai, Bi, Ci = a;Y + bi Y reg + Ci, [= 1,...,6. (170) 


reg 

The equations were used to forecast process air outlet temperature Tout and absolute 
humidity Yout as a function of the temperature and humidity ratio of the regeneration air 
(Tieg Yreg) and of the inlet state of process air (Tin Yin). There are 54 parameters (a; b; c;) 
that related to outlet condition Tout and Yout in the model, so the model is named as 
“Model 54”. Model 54 was carried out to predict the performance of silica gel desiccant 
wheel (type -I). It was observed that the calculated values obtained using this model are 
very close to the measured values. 
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However, for a different type of desiccant material, a new group of 54 coefficients are 
needed, which is a tedious job for users. Then, further simplification were attempted to 
overcome the complexities of “Model 54”. The new model was expressed by the 
parameters ® (relative humidity) and H (enthalpy) 


A® = (Pin — Pour) = A Pin — Prog) +b, (171) 


AH = (Hin — Hou) = a'(H reg — Hin) + 8, (172) 
where the enthalpy is calculated by 


_ 2301+ 1.85T)Y 
7 1000 

the calculation of relative humidity ratio is based on empirical relation developed by 
Beccali [65] 


® = (18.6715 Y + 1.7976)e 037, (174) 


H 


+ 1.006T (173) 


The new relational expressions only have four variables a, b, a’, b’, much less than the 
variable numbers of “Model 54’. Three different kinds of desiccant wheels were used to 
validate this model, two are based on silica gel and another one is on lithium chloride. The 
results showed that the simplified model reasonably predict the performance of all wheels. 

The above two empirical models are valid only for the wheel running with identical air 
flux in process and regeneration side. By adding correction factors, Beccali [66] further 
modified the model which can be used to forecast the system with unequal process and 
regeneration air-flow volume. The modified model was verified for the three different kinds 
of desiccant wheels. 


6. Conclusion 


Various mathematical models have been proposed to assess the performance of 
desiccant wheel for different operating design conditions. Generally, the mathematical 
models are based on heat and mass transfer phenomena inside the system. They can be 
divided into two main categories in terms of gas-side resistance (GSR) model and gas and 
solid-side resistance (GSSR) model. Within the two categories, the models can be 
subdivided based on the dimensions and the methods to represent solid-side resistance 
(SSR). GSR model has two dimension types: one-dimensional and two-dimensional. 
GSSR model can be classified into three sub-categories: pseudo gas-side (PGS) model, gas 
and solid-side (GSS) model and parabolic concentration profile (PCP) model. Also, using 
the measurements data, several empirical models to evaluate the performance of desiccant 
wheel are reviewed in this paper. 

In the GSR model, heat and mass transfer within solid desiccant is not taken into 
account. The governing equations have relatively simple form with lower accuracy. PGS 
model adopts the lumped heat and mass transfer coefficients to represent the overall heat 
and mass transfer from the air stream to desiccant material. In theses models, convective 
heat and mass transfer coefficients are empirically degraded to account for both gas-side 
and solid-side resistances. Therefore, though its form is basically the same as GSR model, 
the accuracy is improved. The more practical model known as GSS model adds second- 
order heat and mass transfer diffusional terms directly into the governing equations. It is 
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more related to the actual process occurred in the desiccant wheel and the effect of SSRs 
can be well explained. The precision of this model is thus greatly improved. PGS model 
consumes less computational time but has lower flexibility and accuracy compared to GSS 
model. PCP model which assumes an analytically convection concentration profile within 
the desiccant particle appears to be a good compromise between fundamental principles, 
accuracy and computational speed. It is an excellent tool for performing seasonal 
simulations or repetitive design simulation particularly when alternative materials are 
investigated. 

It can be seen that boundary conditions in Mode I are widely used in all the models 
whenever possible. Mode II and HI boundary conditions are usually utilized in GSSR 
model. Other auxiliary relations such as heat of sorption and equilibrium isothermal have 
to be obtained to complete the mathematical model. 

Due to the complex nature of the partial differential equations, it is necessary to solve 
the mathematical model using numerical methods. Finite difference method is recognized 
as the most accurate and most universal solution technique. Though other methods such as 
finite volume method [7,38], analogy method [31] and wave analysis method [33,57] are 
also used. 

Beside, model validation is an important step in model development since it offers the 
possibility of comparing simulation results with actual system behavior. Experiments are 
mostly used to validate the mathematical model. Beside, compared with the previous 
numerical results is also a good validation method. 

Though a large amount of works have been conducted on modeling and analyzing 
desiccant wheel, further efforts are still needed: 


(i) To take account of the pressure and velocity loss through the channel. 

(ii) To study the influence of variable thermophysical physical properties of moist air and 
desiccant materials which are taken as constant in almost all existing models. 

(iii) To improve the accuracy of the model based on diffusion mechanism, since by now, 
there still no uniform theory exists. 

(iv) To construct a simulation package of the whole desiccant air-conditioning system by 
combining the model of the wheel and other components in the system. Using object- 
oriented languages, Chau and Worek [67] simulated an open-cycle desiccant cooling 
system operating in the recirculation mode. As presented in the paper, the simulation 
package is ongoing and in the near future it is hoped that a graphical user interface can 
be applied in it (Tables 2 and 3). 
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